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Abstract 

We study the dephasing of a single qubit coupled to a bosonic bath. In particular, we investigate the case when the bath 
is initially prepared in a pure state known as the Schrodinger cat. In clear contradistinction to the time-evolution of an 
initial coherent state, the time evolutions of the purity and the coherence factor now depend on the particular choice of the 
Schrodinger cat state. We also demonstrate that the evolution of the entanglement of a two-qubit system depends on the 
initial conditions in a similar way. 
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1. Introduction 

Controlling the dynamics of open quantum systems 
is of crucial importance for the quantum information 
processing [1]. As there is no general method for ana- 
lyzing the non-Markovian reduced dynamics, the ex- 
actly solvable models may provide important and un- 
biased results. One of the examples is the dephasing 
model [2,3,4,5,6] that describes an idealized case when 
the quantum system does not exchange the energy with 
its environment. This model has recently been studied 
in the context of entanglement dynamics [7,8,9,10] and 
the geometric phases [11]. In particular, it has been 
shown in Ref. [8] that the entanglement can effectively 
be controlled by an external finite bosonic quantum 
system prepared in so-called non-classical states [12]. 

In this paper we study the complementary case when 
the infinite bosonic system is initially prepared in the 
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Schrodinger cat state. For a finite bosonic system such a 
state is defined as a superposition of two coherent states 
with the same amplitudes but with phases shifted by 
TT [12] . Here we generalize this notion to the case of in- 
finite dimensional systems composed of bath and sys- 
tem dynamics. We show that the reduced dynamics 
of the qubit depends on a specific choice of the initial 
Schrodinger cat state. This is in clear contrast to the 
situation when the initial state is purely coherent. It 
holds true not only for purity and coherence of a single 
qubit but also for entanglement of a two-qubit system. 

Due to the decoherence phenomenon, the assumed 
initial state of an infinite bosonic bath is inaccessible in 
the present experiments. However, the development of 
experimental techniques allows one to manipulate and 
control systems devised from an increasing number of 
particles [13]. Therefore, the results presented in this 
paper may serve as a starting point for understanding 
of qubits coupled to large bosonic systems prepared in 
a desired quantum state. Our choice of the initial state 
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is motivated by the fact that multiple Schrodinger cat 
states can accurately approximate any quantum state 
[14,15]. 



2. Model 

We consider a qubit Q, which interacts with the envi- 
ronment R. The Hamiltonian of the total system reads 
[2,3] 



H = Hq®Ir + ] 



i Hr + Hi, 



(1) 



where Iq and Ir are identity operators in correspond- 
ing Hilbert spaces of the qubit Q and the environment 
R, respectively. The qubit Hamiltonian Hq is in the 
form 



Hq =eS^^e{\l}{l\- 1-1} 



(2) 



where the canonical basis of the qubit is {|1), | — 1)} 
and ±e are the energy levels of the qubit. When Q 
represents a particle of spin S = 1/2, the energy e is 
proportional to the magnitude of the external magnetic 
field. The environment is assumed to be a boson field 
described by the Hamiltonian 



Hr = I dui h{Lj)a\Lj)a{Lj), 
Jo 



(3) 



where the real-valued dispersion relation h{uj) speci- 
fies the environment, e.g., h{uj) = u> describes phonon 
or photon environment. The operators a^(Lj) and a(uj) 
are the creation and annihilation boson operators, re- 
spectively. The coupling of the qubit to the environ- 
ment is described by the Hamiltonian 

Hi = \l){l\^H+ + \-l){-l\^H-, (4) 
with 

(5) 



H± = ± dujG{oj)\a{oj) + a\uj)\^ 



where the function G{u!) is the coupling strength. 

Without loosing generality, wo assume that it is a real 
function. The Hamiltonian (1) can be rewritten in the 
form 



H = \l){l\®Hi + \-l){-l\®H.i, 
H±i =Hr + H±±s. 



(6) 
(7) 



Since there is no energy exchange (i.e. we use a non- 
demolition coupling) between the qubit and the en- 
vironment, our modeling corresponds to pure dephas- 
ing. Hamiltonians like (6) have been exploited for de- 
scription of the inter-conversion of electronic and vi- 
brational energy [16], the electron-transfer reactions 
[17], a quantum kicked rotator [18], chaotic dynamics 
of a periodically driven superconducting single elec- 
tron transistor [19] and the Josephson flux qubit [20], 
to mention but a few. 



3. Exact reduced dynamics 

The model we study is exactly solvable [2,4,6], i.e., 
the Schrodinger equation for the wave function 
of the total system can be solved exactly. Here we follow 
the method presented in Ref. [9]. First, one needs to 
specify an initial state |^(0)). Let us assume that at 
the initial time t = 0, the wave function has the form 



|*(0)) = (6i|l>+6_i|-l»®|i?), 



(8) 



where bi and b-i determine the qubit initial state and 
\R} is the initial state of the environment. Then 

Mt)) = bi\l) ® \Mt)) + - 1) ® IV'-iW), (9) 

where \tpi{t)) = exp[—Hit]\R) {i = ±1) can be rewrit- 
ten in the form [9] 

IV'iW) =e-'^^Wz)(s+ e-'««*|i?>, 

=e-'^-i(*)i3(ff- -S,-) e-'"^*\R). (10) 
The phases Ai(t) and A_i(t) are given by 



Ai/_i(t) = ±et 



-f 

Jo 



d(jjg (w) {h(oj)t — sm[h{ijj)t\} , 

(11) 



where the abbreviation g{u>) = G{ui)/h{ui) has been 
introduced. For any function /, the notation ft stands 
for 



-ih{u})t 



(12) 



For an arbitrary square-integrable function /, the dis- 
placement operator D{f) is defined as below [21] 

D{f) = exp dLO [f{oj)a^ioj) - /*(a>)a(a;)] | .(13) 
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The reduced qubit dynamics can be obtained for any 
factorisable initial state of the form 



^(0)= E 



(14) 



where g{0) is the initial statistical operator of the total 
system and pij are non-negative parameters. The re- 
duced statistical operator p{t) for the qubit alone can 
be obtained by tracing the environment degrees of free- 
dom, namely, 

p{t)=TrR [git)] 

J2 PiAi){j\ ® T^H {e-'"**\R}{R\e'"^*) 



i,j=l,-l 

= PijCjiit)\i)(j\, 

i,3=l,-l 



(15) 



where Ttr denotes the partial tracing over the envi- 
ronment variables. Hi for i = ±1 is given by Eq. (7) 
and Cji(t) = {V'j(i)IV'i(i)) is a scalar product between 
the functions \tpj{t)) and |V'i(^)) in the environmental 
Hilbcrt space. The initial state of the qubit \6,<j)} is 
commonly parametrized by two angles on the Bloch 
sphere: The polar angle 6 and azimuthal angle 4>. Then 



\e, 4>) = cos(6l/2)|l) + c^'* sin(6l/2)| - 1). 



(16) 



In this parametrization hi = cos(^/2) and b-i = 
e^'* sin(0/2) (see Eq. (8)) and the initial density matrix 
p(0) of the reduced qubit dynamics reads 



p(0) = 



cos^(6l/2) (l/2)sin6le-*'*' 
sin6le*'* sin^(6l/2) 



\^(l/2) 



(17) 



From Eq. (15) we obtain the density matrix p{t) in the 
form 

/ cos^(0/2) (l/2)A(t)smee^"''' 

Pit) =\ ./ V I -(18) 

\il/2)A* it) sinee"'' sin-(0/2) 

All information about influence of the environment on 
the qubit is incorporated in the dephasing function 
Ait) = c_i,i(t). 

In the following we assume that initially the envi- 
ronment is in the pure Schrodinger cat state, which is 
defined by the relation 



(19) 



where |a) = £'(a)|f2) is the coherent state determined 

by the function a = a(a;) and is the vacuum state 
of the bosonic bath. The normalization constant 



iV = 2 -I- 2cos($)exp 



-2 / doj\aiw)f . (20) 
Jo J 



The phase <I> allows to manipulate the initial state of 
the environment. In this case, the dephasing function 
becomes 

Ait) = iV-^ [(a-i(t)lai(t)) + {Q_i(t)| ~ Qi(t))e** 



+ (-a_i(i)|ai(i))e-'* + (-a_i(t)| - ai(t))J 

(21) 

with |Q±i(t)) = exp(— iJ^±ii)|Q). For the sake of 
brevity we calculate the explicit form of the dephas- 
ing function Ait) for the case of given coherent states 
\a) determined by real functions q(c<;) only. As a first 
main result we find 

Ait) = Af-'Ao(i)e-^'^* |^+(t)e-'* + A_(t)e'* 

-h2cos[4A„(t)]}, (22) 

where 

A„(t)=/ dLuaiiu)gioj)smihioj)t), (23) 
Jo 

ylo(t) =exp|-4^ da;/(a;)[l-cos(/i(a;)t)]|, (24) 
A±(t)=exp|— 2 J dijjQ^i(jj) 

=F4^ dwa(a;)ff(a;)[l-cos(7i(w)t)]|. (25) 

As we show next, the dephasing function Ait) deter- 
mines certain quantifiers describing various aspects of 
quantum information. 

4. Purity and coherence 

We start with basic quantifiers describing the infor- 
mation loss of the qubit. The first one is the purity 
defined by: 

Vit) = Tri/ii)) = \ {\Ait)\^ - 1) sin^ e + 1. (26) 

Its interpretation is clear: The environment results in 
a decrease of the purity. It is equal to 1 for pure states 
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and 1/2 for maximally mixed states. To quantify co- 
herence, we introduce the coherence factor C{t) which 
is determined by the evolution of the non-diagonal el- 
ements of the qubit reduced density matrix, 

|pi2W|=C(i)|pi2(0)|. (27) 
Comparison of (17) and (18) yields 

C{t) = \A{t)\. (28) 

The coherence factor is maximal in the absence of the 
qubit-bath interaction, i.e., C{t) = 1, and vanishes for 
the case of complete decoherence, C{t) = . 

5. Entanglement decay 

In order to study the influence of the dephasing 
on the quantum non-locality, we extend the previous 
model and include a second, completely independent, 
qubit q. The Hamiltonian of such a composite system 
thus reads: 

H=[HQ+HR + Hi]®I^ + Hg, (29) 
Hg = lQiS)lRiS,eS^. (30) 

We assume that the correlations between both the 
qubits are encoded in their initial entanglement. For 
simplicity, we take the depolarized Bell states as the 
initial state, i.e., 

p(0) = (l-p)pi + |lQ®I„ i = l,...,4 (31) 
with 

Pi/2 = ^ [| - 1, 1) ± |1, -1)] K-1, 1| ± (1, , (32) 

P3/4 = ^ [I - 1, -1) ± 11, 1>] K-1- -1| ± ih 1|] . (33) 

The depolarization accounts for an imperfect preparar 
tion of the initial state. 

In a general case, the state of an open system is 
mixed. To quantify its entanglement, several useful 
measures have been proposed [22] . One of the operar 
tional measures is the negativity, defined by N{p) = 
max(0, Ai) , where A; are the negative eigenvalues 
of the partially transposed density matrix of two qubits 
[23]. For the model under consideration the negativ- 
ity can straightforwardly be evaluated for an arbitrary 
evolution time t. One obtains 



JV(p(J)) = max(^0,^|^(t)i-|J. (34) 

The negativity is positive for an entangled mixed state, 
whereas it vanishes for unentangled states. Moreover, 
it presents an entanglement monotone and can be used 
to quantify the degree of entanglement. 



6. Discussion 

The main quantifiers like the purity (26) , the coher- 
ence factor (28) or the negativity (34) depend directly 
on the dephasing function A{t). Therefore, we start 
with discussing its properties in further detail. The 
dephasing function depends on the qubit -environment 
coupling via the functions g{uj) = G{u!)/h{uj) and 
a{u!). The latter one defines the initial Schrodinger 
cat state. For convenience, we can assume that both 
functions are real. We also introduce the new function 
J{u!) = uPg^{(jj). Then, the comparison of the function 
Ao{t) (see Eq. (24)) with the standard expression for 
the decoherence function (see e.g. Eq. (4.51) in Ref 
[6]), allows one to identify J{oj) as the spectral density. 
In the literature, there are several examples of J{io) 
in use. A frequently used one is the generalized Drude 
form defined by [2] 

J{lo) — A oj"'"^'^ exp(— oj/ojc), (35) 

where /x > — 1 and Wc is the cut-oiT frequency. The case 
fj, e (—1, 0) corresponds to a sub-Ohmic, n = Q to the 
conventional Ohmic and ^ € (0, oo) to a super-Ohmic 
environment. 

One can observe that the long-time limit is given by 

Ao = hm Ao{t) = exp |-4 j dwJ(w)/a;^| . (36) 

The integral in this expression is infinite for a sub- 
Ohmic and an Ohmic environment. Then Aq = 
and the dephasing function diminishes to zero, 
linit^so A{t) = 0. Consequently, the purity (26), co- 
herence factor (28) and negativity (34) asymptotically 
take on the following asymptotic long-time values 

p = 1 - 1 sin^ e, C = 0, N = 0. (37) 

One can see that for the sub-Ohmic and Ohmic en- 
vironments all the quantifiers are independent of any 
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particular choice of \a). It means that in the long-time 

regime the qubit properties do not depend any longer 
on the initial Schrodinger cat state. The super-Ohmic 
case is more intriguing because Ao > 0. As it follows 
from the expression for the normalization constant N 
(see Eq. (20)), the function a{u!) is square-integrable. 
Starting from the Cauchy-Schwarz inequality one can 
find that also the integrals in (23) and (25) exist, are 
finite and their values depend on the function a{Lo). In 
consequence, the dephasing function depends on both 
a and $, i.e. on the initial state of the environment. 
Therefore all characteristics (26), (28) and (34) do de- 
pend on the initial environment state, provided the en- 
vironment is super-Ohmic. 

It is instructive to compare Eq. (22) with the de- 
phasing function obtained for an initial, purely coher- 
ent state \R} = I a) of the environment. In this case the 
dephasing function A{t) becomes 

A{t) = exp {-2iet) exp[-4iA„(t)]Ao(t). (38) 

In clear contrast to the initial Schrodinger cat state, 
|^(f)| now depends only on Ao(t). This fact implies 
that the purity (26), the coherence factor (28) and the 
negativity (34) are independent of the initial state of 
the environment also for a super-Ohmic bath. 

7. Conclusions 

Dephasing characteristics of qubits coupled to a 
bosonic environment and prepared in a Schrodinger 
cat state has been investigated. The properties of the 
reduced dynamics, as reflected in the purity and co- 
herence factor, have been shown to exhibit an explicit 
phase-dependence 3> as a parameter of the Schrodinger 
cat state. Qualitatively the same behavior has been 
obtained for the entanglement feature, being quan- 
tified by the negativity. The main conclusion is the 
following: If the initial state of the environment is the 
coherent (or vacuum) state then the informational 
quantifiers do not depend on the initial state. How- 
ever, if the initial state is a linear combination of two 
coherent states as the Schrodinger cat state then such 
quantifiers as the purity, coherence factor or the negar- 
tivity do depend on the initial state, - via the function 
a and the phase $ -, of the environment, at least in 
the short-to-intermediate time regime. For the super- 
Ohmic environment this result holds true also in the 



long-time limit. Moreover, the ^-dependence allows 

one to selectively control the dephasing characteristics 
and the entanglement characteristics. 
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